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Abstract. - Numerically exact path-integral Monte Carlo data are presented for TV < 10 
strongly interacting electrons confined in a 2D parabolic quantum dot, including a defect to 
break rotational symmetry. Low densities are studied, where an incipient Wigner molecule 
forms. A single impurity is found to cause drastic effects: (1) The standard shell-filling sequence 
with magic numbers TV = 4,6,9, corresponding to peaks in the addition energy A(N), is 
destroyed, with a new peak at TV = 8, (2) spin gaps decrease, (3) for TV = 8, sub-Hund's rule 
spin S — is induced, and (4) spatial ordering of the electrons becomes rather sensitive to 
spin. We also comment on the recently observed bunching phenomenon. 



During the past few years, the electronic properties of 2D quantum dots have been the 
subject of intense theoretical studies QUI]. ^ particularly interesting aspect of such artificial 
atoms is the wide experimental tunability of the electron number TV and the Brueckner inter- 
action strength parameter r s in high-quality semiconductor dots 3,4,6 . In particular, it has 
been established both theoretically |1I7| and experimentally 001] that a "Wigner molecule" of 
crystallized electrons forms already at surprisingly high densities corresponding to r s > 2 . 
In this paper, we present numerically exact path-integral Monte Carlo (PIMC) JI] simulation 
results for TV < 10 electrons in the most challenging incipient Wigner molecule regime, choos- 
ing r s ss 4. Due to the breakdown of effective single-particle descriptions, this parameter 
regime is notoriously difficult to treat within approximation schemes such as Hartree-Fock 
theory pQ, density functional theory [Tl ll2l[T3*| . semiclassical or classical ^5] approaches. 
Especially concerning disordered dots, in the absence of benchmark calculations, their accu- 
racy has so far been largely unclear. 

The incipient Wigner molecule regime also exhibits unexpected and interesting physics, in 
particular rather dramatic effects when including just one impurity in order to break rotational 
symmetry. Below we shall demonstrate this in three different ways. First, despite the presence 
of strong interactions, the clean system still exhibits shell structure, with exceptional stability 
of the dot for "magic numbers" TV = 4, 6, and 9. The stability is reflected in peaks in the 
corresponding addition energy, 



A(TV) = E(N + 1) - 2E(N) + E(N - 1), (1) 
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Table I - Quantum dot energies for fixed N and S at UbT = 0.125 from PIMC simulations with and 
without impurity. 
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where E(N) is the energy of the TV-electron dot 0- Strikingly, a single impurity is able to 
completely alter this shell structure. Only the filled-shell peak at TV = 6 survives, and a new 
peak at TV = 8 emerges. A similar peak has also been observed in recent experiments A 
second example concerns the ground-state spin S, where a high value of S for certain TV can 
lead to a "spin blockade" for transport through the dot ^J. The role of impurities on the 
spin polarization of interacting quantum dots has recently been studied by various groups, 
suggesting either a tendency towards spin polarization |T7[. depolarization [TQ, or strongly TV- 
dependent behaviors without general trend |13) . As energy differences between different spin 
states are typically very tiny, approximations are particularly prone to predict false spin. We 
find that "spin gaps" between the ground state and the next-higher spin state are typically 
decreased by the impurity, indicating a trend towards polarization, but without change in 
the ground-state spin for the iV under study. However, there is a notable exception to this 
rule: For TV — 8, the rare case of sub-Hund's rule spin S = is found in the presence of the 
impurity, while the clean dot has S = 1. Finally, as third example, we show that the electron 
crystallization pattern shows a strong dependence on total spin S, especially if rotational 
symmetry is broken. To give a typical example, for TV — 6 electrons with S = 0, electrons 
basically arrange on a ring, but for the fully polarized S = 3 state, the amount of charge in 
the central region increases by about 40% compared to S = 0. Therefore, there is no clear 
separation of energy scales for spin and charge fluctuations. Such behavior is only found in 
the incipient Wigner molecule regime, and invalidates the commonly used lattice-spin models 
appropriate in the deep Wigner-crystallized limit |14| . 

We start by discussing the model underlying our work. A closed 2D quantum dot with 
parabolic confinement of frequency ujq and an additional impurity potential has been studied. 
Measuring energy in units of haj Q and distance in oscillator lengths l a = y/h/m*ijjQ, the 
Hamiltonian for the TV-electron system at zero magnetic field is 

7 = 1 V ' i<j 1 31 



B. Reusch and R. Egger: Impurity effects in few-electron quantum dots 



3 



- a 

X 

'. \ 

2.5- \\ 



1.5- 



— h- 



s s f 



o - 



8 8 



5 

N 



0.2 
0.1 




Fig. 1 - PIMC results with impurity (squares) and for a clean dot (circles), for a) the addition energy 
A(iV) and b) the spin gap 8E(N) for N < 8. Error bars, unless shown explicitly, are of the order of 
the symbol size. Dashed and dotted lines are guides to the eye only. 



where A = lo/a for effective Bohr radius a. Throughout this study, we set A = 4, which puts 
the quantum dot into the incipient Wigner molecule regime with r s « 4 pQ. Here r s is half 
the nearest-neighbor distance of the electrons measured in units of a. For GaAs-based dots, 
this corresponds to Hluq pa 0.7 meV 015] • We allow for an attractive impurity of strength w 
centered at x. The impurity potential is taken in the form |13j 

V imp (r) = -u,exp[-2(r - xf/a 2 ], (3) 

with parameters chosen as w — 4, x — (0, 1.5), and impurity range a = 0.75. Qualitatively 
similar results were also observed for other parameter choices. For the above parameters, the 
impurity is close enough to the center to influence the electrons. Typically, when integrated 
over a circle of radius a around x, about 3/4 of an electron charge is trapped by this impurity. 
However, the impurity location is also sufficiently far away from the center to break rotational 
symmetry. 

Up to statistical error bars, PIMC data presented below are numerically exact. For A = 4 
and fcflT = 0.125, the fermion sign problem 1 1 1 j permits us to extract reliable results only for 
N < 10. The simulations yield accurate estimates, at given N and total spin S = (iVf — iVjJ/2, 
for the energy, the charge or spin density, and other correlation functions of interest. Typically, 
5 x 10 7 MC samples were accumulated for each parameter set. Errors due to finite Trotter time 
discretization At were carefully eliminated by extrapolating results for different discretization 
down to At — > with a linear regression fit, using the fact that the expectation value of any 
Hermitian observable has finite- At corrections vanishing tx Ar 2 |19| . This extrapolation 
procedure is essential for achieving high accuracy, as Trotter convergence is very slow except 
for fully spin-polarized clean dots [231 ■ PIMC results for the spin-dependent energy are listed 
in table[H where bracketed numbers denote one standard deviation error bar for the last digit, 
which includes stochastic MC errors and extrapolation uncertainties. When comparing our 
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Fig. 2 - Shadowed contour plots of the density (with impurity) for N = 7,8,9, and 10 [a)-d)]. 
Different contours lie at integral multiples of 0.1 times the maximal density; darker colors signal 
higher density. In this figure, ksT = 0.25, with very similar but noisier results at lower T. 



data in table [I] to T = results obtained from other methods such as exact diagonalization 
|21I22I23") . one has to take into account finite temperature, which can cause sizeable differences 
for not fully spin-polarized states. For instance, for two electrons and S = 0, PIMC predicts 
Mclean = 4.893(7), while the T = result is E = 4.848. Such differences arise because 
PIMC proceeds at fixed total z-component of the spin and therefore receives contributions 
from higher spin states [Jj. Note that this effect is rather tiny and does neither affect the 
(finite-temperature) addition energies Q nor the prediction of ground-state spins. 

Let us first discuss the addition energies displayed in fig. [TJi). Remarkably, in a clean dot, 
one can still observe the standard shell-filling sequence also found experimentally in small 
vertical dots |4ll2*4]. The filled-shell configuration N = 6 and the half-filled cases N — 4, 9 are 
exceptionally stable. The peaks in the addition energy at these magic numbers signal that the 
Wigner molecule is not yet fully established, since the classical Wigner limit is characterized 
by different and much less pronounced magic numbers [T^. By simultaneously monitoring the 
charge density (see below), however, clear traces of spatial ordering can be observed, reflecting 
the onset of Wigner crystallization. The energetic shell structure is then destroyed once the 
impurity is present, see fig.^Ji). Apart from a new peak at N — 8, the addition energies A(iV) 
now show a monotonic decrease with N. Only the filled-shell peak at N = 6 is still visible, yet 
considerably diminished. We note that the peak at N = 2 is always absent for this interaction 
strength Q]E|. 

Next we address spin properties in this regime, see tableland fig-CJ)). Except for N = 8, 
the spin S was not altered by the impurity |25|. While S mostly takes the minimal value 
(5 = or S = 1/2), for N = 4 {N = 9), we have Hund's rule spin S = 1 (5 = 3/2). 
Furthermore, for N — 10, also a partially spin-polarized state S = 1 is found. For the clean 
N = 3 dot, at A c w 4.34 a spin transition S = 1/2 -> 3/2 occurs at T = 0II1I221- We 
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Fig. 3 - Contour plots of the density for N — 6. a) With impurity, S = 0. b) Clean case, S = 0. c) 
With impurity, S — 3. d) Clean case, S = 3. 



are already very close to this transition point, and within error bars, we cannot decide the 
ground state spin. The case N = 8 is particularly noteworthy. For the clean dot, the spin 
is S = 1 according to Hund's rule, but the S = 2 state is actually very close, see tabic [I] 
Strikingly, with impurity, we then find sub-Hund's rule spin S — 0. Since this effect requires 
sufficiently strong interactions (low density), for a weakly disordered dot with N — 8, one 
can expect a transition from the usual Hund's rule S = 1 ground state at high density (small 
r s ) to the unpolarized S — state at low density. Since for TV = 9, we have spin S = 3/2, 
this transition could be observed experimentally as spin blockade for the 9th electron to enter 
the dot. By monitoring spin-spin correlations between the electrons at the impurity location 
and the central region (data not shown here), strong antiferromagnctic couplings are visible. 
Hence it is likely that local moment formation at the impurity location and in the central 
region is important for this sub-Hund's rule behavior. This is also confirmed by noting that 
charge and spin densities do not follow each other, as they do in the deep Wigner-molecule 
regime jUEHl- Since except for N = 8 the spin S is not changed, the impurity's overall effect 
on the spin polarization seems small. Nevertheless, our results show a decreasing "spin gap" 
5E(N) = E(S+1) —E(S), see fig-H 3 ), indicating a weak tendency towards polarization. Note 
that the spin excitation energies are very small, again stressing the need for highly accurate 
calculations. 

We now turn to the onset of electron crystallization in real space. Figure shows density 
plots for N — 7 to 10 electrons in their respective ground-state spin configuration. At the 
defect, the density is about three times larger than elsewhere in the dot. Although there 
are large quantum fluctuations, a spatial shell structure is already discernible. The classical 
prediction for clean dots at these electron numbers is as follows ^Hl- For N = 7 and N = 8, 
just one electron is in the center, while for TV = 9 and N = 10, the spatial filling is (2 — 7) 
and (2 — 8), respectively, with the configurations (1 — 8) and (3 — 7) very close in energy. 



6 



EUROPHYSICS LETTERS 



0.45 
0.4 

o 

c 0.35 
0.3 
0.25 



N=6 



Fig. 4 - Amount n c of electron charge localized in the center (see text) for N = 6 and different S, 
both with (squares) and without (circles) the defect. 



Quantum-mechanically, for A = 4, the radial ordering is not yet very pronounced in the clean 
dot |7ll26| . With impurity, however, our data in fig. [5] displays spatial shells, although with a 
rather different filling sequence. Detailed examination of the central region shows that only 
one to two electrons are accommodated in this part. Therefore incoming electrons enter the 
outer ring for N = 7, 8, 9, and 10, with the ring expanding and, for N = 10, nearly embracing 
the impurity. While for very large r s , the crystal structure is S'-independent and effective 
Heisenberg-type lattice models for the spin physics can be constructed, for r s w 4, the spatial 
charge ordering is strongly dependent on total spin S; see also [531|2Sj. This is seen in fig. 
for N = 6, where the densities for S — and S = 3 are compared, both with and without 
the defect. Interestingly, for 5 = 0, electrons tend to basically arrange on a ring, while for 
the spin-polarized case, the classical configuration with one electron in the center |15U2fi| 
is approached. This effect, in particular concerning the influence of the impurity, can be 
quantified by computing the fraction n c of electron charge localized within a radius a = 0.75 
around the center of the dot. The result is shown in fig. 0] and illustrates that more and 
more charge goes into the central region when S is increased. For the clean dot, comparing 
S = and S = 3, this amounts to a w 25% increase in n c , but with the defect, this increase is 
« 40%. That implies that the spin sensitivity of the incipient Wigner crystallization process 
is enhanced by the impurity, probably due to a local spin at the impurity site. 

Finally, let us briefly discuss the relation of our results to the experimentally observed 
"bunching" phenomenon found in the Wigner molecule regime where the addition 

energy was found to vanish at certain values of N. Bunching was also reported in a recent 
numerical diagonalization study for N — 7 and N — 8 [2E|> where the dot was modeled by a 
3x4 Hubbard- type model containing only short-range interactions. In contrast to these results, 
for the more realistic continuum Hamiltonian (0) with long-range interactions, bunching is 
not observed for N < 10, see fig. QJi). Since bunching is believed to be a generic effect [27] . 
a consistent explanation probably has to assume strong disorder or other mechanisms not 
contained in our model. 

To conclude, we have presented a PIMC study of impurity effects in 2D quantum dots. We 
have focused on the most difficult yet most interesting regime of an incipient Wigner molecule. 
Large effects are caused already by a single impurity. When the rotational symmetry of the 
clean dot is broken, magic numbers for the energetic shell-filling sequence, observed as peaks 
in the addition energy, are changed. Instead of the usual peaks at N = 4, 6, 9, with impurity 
only a (weak) filled-shell peak at N = 6 and a new peak at N = 8 are observed. In addition, 
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for N = 8, we find a sub-Hund's rule spin 5 = 0. Finally, the incipient Wigner crystallization 
process was found to be strongly spin-dependent, in particular when the impurity is present. 
We hope that these results will also be valuable as reference data to assess the validity of 
computationally less demanding yet approximate calculations. 

* * * 
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